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Abstract 

A long flexible neutral polymer chain immersed in a poor solvent and interacting 
with an impenetrable attractive surface exhibits a phase known as surface attached 
globule (SAG) in addition to other adsorbed and desorbed phases. In the ther- 
modynamic limit, the SAG phase has the same free energy per monomer as the 
globular phase, and the transition between them is a surface transition. We have 
investigated the phase diagrams of such a chain in both two- and three- dimen- 
sions and calculated the distribution of monomers in different domains of the phase 
diagram. 
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A long flexible neutral polymer chain immersed in a poor solvent and inter- 
acting with an impenetrable surface is known to exhibit a very rich phase 
diagram [1]. Competition between the lower internal energy near an attrac- 
tive wall and higher entropy away from it results in a transition, where for 
a strongly attractive surface the polymer sticks to the surface and for weak 
attraction the polymer chain remains desorbed. On the other hand, due to the 
self-attraction in the polymers the possibility of a collapse transition both in 
the desorbed and adsorbed states occur. The phase diagram is generally plot- 
ted in variables uj = e^" and u = e^ eu where (3 is the inverse temperature, e s 
is the (attractive) energy associated with each monomer lying on the surface 
and e u represents an attractive interaction energy between pairs of monomer 
which come close to each other and are separated along the chain by more 
than one unit. The short range repulsion between the monomer is taken in a 
lattice model by self-avoidance [2]. Considerable attention has been paid to 
find the full phase diagram of such a chain in both two dimensions (2D) and 
three dimensions (3D). 

In one of the earliest papers on the subject, Bouchaud and Vannimenus [3] 
derived the exact phase diagram on a 3D-Sierpinski gasket. The phase dia- 
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Fig. 1. Schematic diagrams of self-avoiding walk conformation at T=0. Below crit- 
ical value w c i,(a) the number of monomers at the surface are zero and polymer is 
in the desorbed state; (b) the structure shown in Fig. (a) gets adsorbed on the sur- 
face at critical value to c i . Average number of pair of neighbouring monomers on the 
surface < N c > scales as N^ 1 ; (c) < Nc > increases with u>; (d) At u> C 2, < N c > 
is maximum; (e) < N c > decreases sharply with further increase in w. 

gram consisted of the adsorbed expanded (AE), desorbed expanded (DE) and 
desorbed collapsed (DC) phases. Kumar and Singh [4] have reinvestigated the 
phase diagram for 3-D-Sierpinski gasket and showed that for a certain range 
of surface interactions, an additional phase having the feature of globule at- 
tached to the surface exists. Recently Singh et.al. [5,6] using extrapolation of 
exact series expansions calculated the phase diagram for the polymer chain 
in both 2D and 3D (Euclidean) space. These phase diagrams show the exis- 
tence of the surface attached globule (SAG) phase in qualitative agreement 
with the one found earlier for the gasket. It is important to note that in the 
thermodynamic limit, the SAG phase has the same free energy per monomer 
as the DC phase, and the transition between them is a surface transition. 

As shown earlier [1], it is easy to understand the SAG and its boundary in 
the limit T = 0. At T = 0, the polymer configurations may be mapped by 
Hamiltonian walk as shown in Fig 1(a). Its bulk energy is the same as in the 
DC phase and the surface energy is (e u — e s )Lj| _1 + (d— \)e u N/L\\. Minimizing 
the surface energy with respect to L», we obtain 



E SAG = ~( d - l ) e " N + de ^ d ( e « - e s )dN— 
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While for the DC phase, 



E DC 



(d 



l)e u N + de u N^ . 
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Comparing the energies of these phases, we see that SAG phase has lower 
free energy than the DC or the AC phases for < e s < e u . Thus the lower 
and upper boundaries of the SAG phase (lines u cl and uj c2 in Fig. 2) tend to 
cu cl = 1 and w c2 = u for large u. 

In the case of partially directed polymer in 2D we have found [1] exact phase 
diagram of the SAG phase. In this case the polymer has different behaviour 
depending on whether it is near the wall perpendicular to the preferred direc- 
tion (SAG1) or the wall is parallel to the preferred direction (SAG2). The 
phase boundaries of SAG1 and SAG2 have been found by calculating this 
orientation dependent surface energy. 

Since the analytical approach are limited to very few cases such as fractal 
lattices or the directed walks on 2D lattices, we have to resort to numerical 
methods to calculate the phase diagram in other cases. A lattice model using 
extrapolation of exact series expansions (herein after referred to exact enu- 
meration method) has been found to give satisfactory results as it takes into 
account the corrections to scaling. To achieve the same accuracy by the Monte 
Carlo method, a chain of about two order magnitude larger than in the exact 
enumeration method has to be considered [7]. 

In this article we report our results which we have found using exact enumer- 
ation method both for 2D and 3D cases. Let Cn(N s , N u ) be the number of 
SAWs of A" sites having N s monomers on the surface and N u nearest neigh- 
bour monomer pairs. We analyzed the series for Cn(N s , N u ) up to A^ = 30 for 
square lattice and A^ = 20 for the cubic lattice; thus extending the series by 
two more terms in both cases. 

To obtain better estimates of critical points as well as the phase boundaries, 
we extrapolate for large A" using ratio method [8]. Let 



be the partition function. Then, the reduced free energy per monomer can be 



Z n (uj,u)= C N (N s ,N u )u N °u 



(4) 
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Fig. 2. This figure shows the phase diagram of a surface interacting linear polymer 
in 3D and 2D space respectively. The transition lines were obtained using method 
of exact-enumeration. 



written as [8] 



G(u,u) = lim ^- log Z n (uj,u). (5) 



The phase boundaries are then found from the maxima of 9 G q^' u ^ (= dl ^ a ^ 

del V de u 



an d FGfau) (_ 8(N U ) 



The phase diagram thus obtained is shown in Fig. 2 (a) for three dimensions. 
It has five phases which are denoted as desorbed expanded (DE), desorbed 
collapsed (DC), adsorbed expanded (AE), adsorbed collapsed (AC) and 
surface attached globule (SAG). Since in case of two dimensions, surface 
is a line, therefore, we do not find adsorbed collapsed phase [Fig2(b)]. The 
value of u c at tu = 1 is found to be 1.93 ± 0.02 and 1.76 ± 0.02 for 2D 
and 3-D, respectively. The value of u c at u — 1 is found to be 2.05 ± 0.01 
and 1.48 ± 0.02 for 2D and 3D, respectively. The special adsorption line u c 
separates the AE phase from that of DE phase. The u c line meets u c line 
which is a multi-critical point (u c ,u c ) = (1.93,1.46) for 2D and (1.76, 1.38) 
for 3D, respectively. Values of 3D are same as those reported by us earlier [1]. 
The value of critical exponents both at the ordinary and multi-critical points 
are in very good agreement with the known results [?]. 

In order to show how the monomers distribute when the chain is in the dif- 
ferent regimes of the phase diagram, we calculate the value of each term in 
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Fig. 3. Plot of £/£ m as a function of N s and N p . Here £ = C N (N S , N p )lu n °u n p and 
£ m is the maximum value of £ for a given value of u> and u. (a) u = 1.0, cj = 1.0; 
(b) tt = 1.0, w = 2.0; (c) u = 1.0, w = 3.0; (d) u = 1.0, w = 3.5 




Fig. 4. Same as Fig.3, (a) u = 1.5, w = 1.0; (b) u = 1.5, uj = 2.0; (c) u = 1.5, 
c<j = 3.0; (d) u = 1.5, uj = 3.5 

Eq. (4) for given values of u and u and plot the results for walk up to 20 steps 
(3D case only) in Fig. (3-6). Here we have plotted £/£ m for eight values of uj, 
where £ = C^(N S , N u )uj Ns u Nu (£ m being the maximum value of £) as a func- 
tion of N s and A^ M . Since the term which contribute most to the sum of Eq.(4) 
represents the most probable configurations of the chains, we can infer, from 
these figures about the manner in which the monomers are distributed in dif- 
ferent states. For example, Fig.(3(a-d) and 4(a-d)) show how the distribution 
of monomers changes as u is varied at fixed u — 1.0 and 1.5 (corresponding to 
expanded state). It is clear from these figures (3(a) and 4(a)) that N s and N u 
are nearly zero around u — 1. In Fig.(3(d)-4(d)), at uj = 3.5, iV s is around 20 
corresponding to the adsorbed expanded state of chain. For u = 2.0 and 2.5 
(correspoding to the collapsed state), distribution is shown in Fig.(5(a-d) and 
6(a-d). For u — 1, the polymer chain is in the bulk as there is any monomer on 
the surface (iV s ~ 0) but N u is large which indicate that monomers occupy the 
neighbouring sites. When uj is chosen equal to 2.0 a value which lies in between 
uj c I and uj c 2 line, as shown in Fig. 2, a fraction of monomers gets attached to 
the surface (iV s ~ 12) indicating the breaking of translational invariance of the 
globule corresponding to the SAG phase. With further increase in uj polymer 
gets adsorbed in collapsed state where N s ~ 20 and N u is around 8. 

We also calculate number of neighbouring pairs of monomers at the surface 
i.e. < N c >. This can be calculated from the expression 

Y\N c Zm{uj,u) , \ 

< n c >= y v 6 

Y,Z n (uj,u) 
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Fig. 5. Same as FIG.3, (a) u = 2.0, u = 1.0; (b) u = 2.0, to = 2.0; (c) u = 2.0, 
u = 3.0; (d) u = 2.0, u = 3.5 




Fig. 6. Same as FIG.3, (a) u = 2.5, w = 1.0; (b) u = 2.5, w = 2.0; (c) u = 2.5, 
w = 3.0; (d) u = 2.5, w = 3.5 
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Fig. 7. Diagram shows the variation of number of monomer pairs on the surface 
< N c > with to. 

The variation of < iV c > with u for different value of u — 10, 15, 20, 25 is 
shown in Fig. 7 for two dimensions. At such a high value of u the most likely 
structure is of Hamiltonian walk as shown in Fig. 1(a). Below certain value of 
u c i, this structure remains in bulk and < N c > and < N s > are almost zero. At 
certain critical value u> c i, this structure gets stuck on the surface (Fig. 1(b)). 
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At this value of uj, the < N c > is equal to < N s > and goes as jV d_1 . As 
uj increases < N c > also increases and reaches to its maximum at uj C 2- The 
corresponding configuration is shown in Fig. 1(d). With further increase in uj, 
polymer gets adsorbed on the surface with decrease in < N c > as shown in 
Fig. 7 and Fig. 1(e). 

The results given above should provide ample evidences in favour of the exis- 
tence of the SAG phase. The phase boundaries lo c \ and uj C 2 tend to 1 and u 
respectively as T goes to zero and is independent of dimension. We once again 
emphasized that the transition associated with uj c1 line is surface transition. 

We thank Deepak Dhar and R. Rajesh for many helpful discussions. We would 
also like to thank the organizers for giving us opportunity to present this 
work. Financial assistance from INS A, New Delhi and DST, New Delhi are 
acknowledged. 
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Abstract 

A long flexible neutral polymer chain immersed in a poor solvent and interacting 
with an impenetrable attractive surface exhibits a phase known as surface attached 
globule (SAG) in addition to other adsorbed and desorbed phases. In the ther- 
modynamic limit, the SAG phase has the same free energy per monomer as the 
globular phase, and the transition between them is a surface transition. We have 
investigated the phase diagrams of such a chain in both two- and three- dimen- 
sions and calculated the distribution of monomers in different domains of the phase 
diagram. 
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A long flexible neutral polymer chain immersed in a poor solvent and inter- 
acting with an impenetrable surface is known to exhibit a very rich phase 
diagram [1]. Competition between the lower internal energy near an attrac- 
tive wall and higher entropy away from it results in a transition, where for 
a strongly attractive surface the polymer sticks to the surface and for weak 
attraction the polymer chain remains desorbed. On the other hand, due to the 
self-attraction in the polymers the possibility of a collapse transition both in 
the desorbed and adsorbed states occur. The phase diagram is generally plot- 
ted in variables to = e@ ts and u = e^ 6 " where j3 is the inverse temperature, e s 
is the (attractive) energy associated with each monomer lying on the surface 
and e u represents an attractive interaction energy between pairs of monomer 
which come close to each other and are separated along the chain by more 
than one unit. The short range repulsion between the monomer is taken in a 
lattice model by self-avoidance [2]. Considerable attention has been paid to 
find the full phase diagram of such a chain in both two dimensions (2D) and 
three dimensions (3-D). 

In one of the earliest papers on the subject, Bouchaud and Vannimenus [3] 
derived the exact phase diagram on a 3D-Sierpinski gasket. The phase dia- 
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Fig. 1. Schematic diagrams of self- avoiding walk conformation at T=0. Below critical 
value w c i,(a) the number of monomers at the surface are zero and polymer is in the 
desorbed state; (b) the structure shown in Fig. (a) gets adsorbed on the surface at 
critical value co c \. Average number of pair of neighbouring monomers on the surface 
< N c > scales as N d ~ lm , (c) < Nc > increases with ui\ (d) At w C 2, < N c > is 
maximum; (e) < N c > decreases sharply with further increase in lj. 

gram consisted of the adsorbed expanded (AE), desorbed expanded (DE) and 
desorbed collapsed (DC) phases. Kumar and Singh [4] have reinvestigated the 
phase diagram for 3-D-Sierpinski gasket and showed that for a certain range 
of surface interactions, an additional phase having the feature of globule at- 
tached to the surface exists. Recently Singh et.al. [5,6] using extrapolation of 
exact series expansions calculated the phase diagram for the polymer chain 
in both 2D and 3D (Euclidean) space. These phase diagrams show the exis- 
tence of the surface attached globule (SAG) phase in qualitative agreement 
with the one found earlier for the gasket. It is important to note that in the 
thermodynamic limit, the SAG phase has the same free energy per monomer 
as the DC phase, and the transition between them is a surface transition. 

As shown earlier [1], it is easy to understand the SAG and its boundary in 
the limit T = 0. At T = 0, the polymer configurations may be mapped by 
Hamiltonian walk as shown in Fig 1(a). Its bulk energy is the same as in the 
DC phase and the surface energy is (e u — e^Ly -1 + (d— l)e u N/L\\. Minimizing 
the surface energy with respect to Ln, we obtain 



E SAG = -(d- l)e u N + de u d (e u - e s )*N 



d-l 

d 



(1) 



2 



While for the DC phase, 



E DC 



(d- l)e u N + deuN^cT . 



(2) 



In the AC 

T = is 



phase, we have L± 



1, L\\ = N l K d l \ and the free energy at 



E AC 



(d-2)e u N-e s N + {d-l)e u N^. 



(3) 



Comparing the energies of these phases, we see that SAG phase has lower 
free energy than the DC or the AC phases for < e s < e u . Thus the lower 
and upper boundaries of the SAG phase (lines u c i and w c2 in Fig. 2) tend to 
uj c i = 1 and w c2 = u for large u. 

In the case of partially directed polymer in 2D we have found [1] exact phase 
diagram of the SAG phase. In this case the polymer has different behaviour 
depending on whether it is near the wall perpendicular to the preferred direc- 
tion (SAG1) or the wall is parallel to the preferred direction (SAG2). The 
phase boundaries of SAG1 and SAG2 have been found by calculating this 
orientation dependent surface energy. 

Since the analytical approach are limited to very few cases such as fractal 
lattices or the directed walks on 2D lattices, we have to resort to numerical 
methods to calculate the phase diagram in other cases. A lattice model using 
extrapolation of exact series expansions (herein after referred to exact enu- 
meration method) has been found to give satisfactory results as it takes into 
account the corrections to scaling. To achieve the same accuracy by the Monte 
Carlo method, a chain of about two order magnitude larger than in the exact 
enumeration method has to be considered [7]. 

In this article we report our results which we have found using exact enumer- 
ation method both for 2D and ZD cases. Let Cn{N s , N u ) be the number of 
SAWs of N sites having N s monomers on the surface and N u nearest neigh- 
bour monomer pairs. We analyzed the series for Cn(N s , N u ) up to N = 30 for 
square lattice and N = 20 for the cubic lattice; thus extending the series by 
two more terms in both cases. 

To obtain better estimates of critical points as well as the phase boundaries, 
we extrapolate for large N using ratio method [8]. Let 



be the partition function. Then, the reduced free energy per monomer can be 



Z N (u,u)= C N (N s ,N u )u N °u 



N, 



(4) 
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Fig. 2. This figure shows the phase diagram of a surface interacting linear polymer 
in 3D and 2D space respectively. The transition lines were obtained using method 
of exact-enumeration. 

written as [8] 

G(u),u) = lim — \ogZ N (u, u). (5) 



The phase boundaries are then found from the maxima of 9 G ^' u>> (= at ^^ ) 



del ^ de 3 

)_\ 

del ^ dt u 



, d*G(uj,u) ,_ d(N u ) s 

ancl ~~a^ I" ~lhT)- 



The phase diagram thus obtained is shown in Fig. 2(a) for three dimensions. 
It has five phases which are denoted as desorbed expanded (DE), desorbed 
collapsed (DC), adsorbed expanded (AE), adsorbed collapsed (AC) and 
surface attached globule (SAG). Since in case of two dimensions, surface 
is a line, therefore, we do not find adsorbed collapsed phase [Fig2(b)]. The 
value of u c at uj = 1 is found to be 1.93 ± 0.02 and 1.76 ± 0.02 for 2D 
and 3D, respectively. The value of uj c at u = 1 is found to be 2.05 ± 0.01 
and 1.48 ± 0.02 for 2D and 3D, respectively. The special adsorption line uj c 
separates the AE phase from that of DE phase. The u c line meets u c line 
which is a multi-critical point (u c ,u c ) = (1.93,1.46) for 2D and (1.76, 1.38) 
for 3-D, respectively. Values of 3D are same as those reported by us earlier [1]. 
The value of critical exponents both at the ordinary and multi-critical points 
are in very good agreement with the known results [?]. 

In order to show how the monomers distribute when the chain is in the dif- 
ferent regimes of the phase diagram, we calculate the value of each term in 
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Fig. 3. Plot of £/£ m as a function of N s and N p . Here £ = C N {N S , N p )oj n "u n p and 
£ m is the maximum value of £ for a given value of co and u. (a) i/ = 1.0, lj = 1.0; 
(b) it = 1.0, w = 2.0; (c) u = 1.0, lj = 3.0; (d) u = 1.0, a; = 3.5 




Fig. 4. Same as Fig.3, (a) u = 1.5, u = 1.0; (b) it = 1.5, u = 2.0; (c) u = 1.5, 
u> = 3.0; (d) it = 1.5, oj = 3.5 

Eq.(4) for given values of u and u and plot the results for walk up to 20 steps 
(3D case only) in Fig. (3-6). Here we have plotted £/£ m for eight values of u, 
where £ = Cn(N s , N u )u Ns u Nu (£ m being the maximum value of £) as a func- 
tion of an d A^ M . Since the term which contribute most to the sum of Eq.(4) 
represents the most probable configurations of the chains, we can infer, from 
these figures about the manner in which the monomers are distributed in dif- 
ferent states. For example, Fig.(3(a-d) and 4(a-d)) show how the distribution 
of monomers changes as u is varied at fixed u = 1.0 and 1.5 (corresponding to 
expanded state). It is clear from these figures (3(a) and 4(a)) that N s and N u 
are nearly zero around to = 1. In Fig.(3(d)-4(d)), at to = 3.5, N s is around 20 
corresponding to the adsorbed expanded state of chain. For u = 2.0 and 2.5 
(correspoding to the collapsed state), distribution is shown in Fig.(5(a-d) and 
6(a-d). For to = 1, the polymer chain is in the bulk as there is any monomer on 
the surface (N s ~ 0) but N u is large which indicate that monomers occupy the 
neighbouring sites. When lo is chosen equal to 2.0 a value which lies in between 
lo c 1 and uj c 2 line, as shown in Fig. 2, a fraction of monomers gets attached to 
the surface (N s ~ 12) indicating the breaking of translational invariance of the 
globule corresponding to the SAG phase. With further increase in u polymer 
gets adsorbed in collapsed state where N s ~ 20 and N u is around 8. 

We also calculate number of neighbouring pairs of monomers at the surface 
i.e. < N c >. This can be calculated from the expression 

< N c >= ^ff-f (6) 
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Fig. 5. Same as FIG. 3, (a) u = 2.0, u = 1.0; (b) u = 2.0, u = 2.0; (c) u = 2.0, 
u = 3.0; (d) u = 2.0, u = 3.5 




Fig. 6. Same as FIG. 3, (a) u = 2.5, u = 1.0; (b) « = 2.5, w = 2.0; (c) u = 2.5, 
u = 3.0; (d) u = 2.5, w = 3.5 
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Fig. 7. Diagram shows the variation of number of monomer pairs on the surface 
< N c > with lj. 

The variation of < iV c > with u for different value of u = 10, 15, 20, 25 is 
shown in Fig. 7 for two dimensions. At such a high value of u the most likely 
structure is of Hamiltonian walk as shown in Fig. 1(a). Below certain value of 
cj c i, this structure remains in bulk and < N c > and < N s > are almost zero. At 
certain critical value u c i, this structure gets stuck on the surface (Fig. 1(b)). 
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At this value of uj, the < N c > is equal to < N s > and goes as N 4 ^ 1 . As 
uj increases < N c > also increases and reaches to its maximum at ui C 2. The 
corresponding configuration is shown in Fig. 1(d). With further increase in uj, 
polymer gets adsorbed on the surface with decrease in < N c > as shown in 
Fig. 7 and Fig. 1(e). 

The results given above should provide ample evidences in favour of the exis- 
tence of the SAG phase. The phase boundaries ui c i and ui C 2 tend to 1 and u 
respectively as T goes to zero and is independent of dimension. We once again 
emphasized that the transition associated with ui c i line is surface transition. 

We thank Deepak Dhar and R. Rajesh for many helpful discussions. We would 
also like to thank the organizers for giving us opportunity to present this 
work. Financial assistance from INS A, New Delhi and DST, New Delhi are 
acknowledged. 
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